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I. INTRODUCTION 

The conventional quantum description of open systems employs density matri- 
ces and their nonunitary evolution equations, i. e. master equations. An equivalent 
description can be given in terms of wave functions, at the expense of introducing 
damping and noise terms in the Schrodinger equation (see e. g. P|-[2^). For a given 



open system, both the master equation and the corresponding stochastic Schrodinger 
equation can be derived by treating the interaction, usually weak, with an environ- 
ment. 

Whether one deals with a master equation or the equivalent stochastic 
Schrodinger equation is mostly a matter of convenience. Less work is involved in 
determining A^ components of a state vector for a single realisation of the noise pro- 
cess rather than A^^ elements of a density matrix; however, many realisations of the 
noise must be followed in order to establish the full information stored in the density 
matrix. 

Apart from occasional computational convenience the stochastic Schrodinger 
equation has a certain conceptional appeal. Recent experiments on individual ions 
or atoms stored in electromagnetic traps have revived the old question about what 
quantum mechanics has to say about the course of a single run of an experiment. 
The so-called quantum jumps of an atom into or out of a metastable state and the 
accompanying random-telegraph modulation of a fluorescence signal, for instance, 
strengthen our interest in a formulation of dissipative quantum dynamics where ran- 
dom jumps are naturally reflected in the time dependence of the wave function. 

The stochastic Schrodinger equation is particularly welcome to some researchers 
interested in quatum measurement theory. The seemingly scandalous effective col- 
lapse of a superposition to a mixture upon measurement is by now rather well un- 
derstood, even on the basis of exactly solvable models p3|-|2^. However, it is still 



esthetically satisfying to see the wavefunction of a pointer variable evolve, in each 
realization of the noise, to one particular distinguished displacement corresponding 
to one eigenvalue of the measured observable of some microsystem; and to see many 
independent realisations of the noise build up the relative frequencies of the possible 
outcomes. 

Other peculiarities of quantum mechanics, like Schrodinger's cat and the EPR 
paradox, appear under an interesting new perspective in the description by a stochas- 
tic Schrodinger equation. The fate of Schrodinger's cat was recently investigated 
numerically in a simple model in this manner p6|. Previous treatments of such 



problems using the density matrix [^,^ were thus complemented in rather more 
intuitive ways. 

We here take up the cat problem again. We present a rigorous solution of a 
stochastic Schrodinger equation for the damped harmonic oscillator, starting from 
an initial state of the cat type. It is most gratifying to see how the wavevector chooses 
randomly between the two possible and "macroscopically" distinguishable locations. 
The choice is made on the so-called decoherence time scale tdec which is shorter than 
the mechanical damping time by a factor measuring the difference between the two 
initial locations. For each realization of the noise, jumps between the two locations 
become exceedingly unlikely after several units of tdec- 

Noise can be accounted for in the stochastic Schrodinger equation in different 
but essentially equivalent ways. One possibility is to let the noise represent vacuum 
and/or thermal fluctuations of the environment; it is then natural to speak of input 
noise. Alternatively, we may let the noisy part of the output of the system appear as 
stochastic driving. Both forms of the noise are of course related to one another, as 



described by the input-output formalism of Collet and Gardiner |^ . Depending on 



which form of the noise is chosen in the stochastic Schrodinger equation, the latter 



takes a linear or nonlinear appearance. The equivalence of the two formulations, 



recently shown by Ghirardi, Pearle and Rimini [|T5|, by Belavkin and Staszewski |T^ 
and by Goetsch and Graham [^,^], will be illustrated again below. 

II. LINEAR EVOLUTION FROM A CAT STATE 



We first discuss the solution of the linear stochastic Schrodinger equation [ll5H2^ 



d\Mt)) = [-^a^adt + ^ae-"^de{t)^ \Mt)) ■ (2.1) 

Here d6{t) is a real noise increment independent of its predecessors in time; otherwise, 
d6{t) need not be specified at this point. The index 6 on the wave vector is meant to 
indicate the contingency of that vector on the noise history: \iljg{t)) is a functional of 
the noise increments arisen in the past but independent of the one arising at present, 
d6{t). To include the description of homodyning with a local oscillator of constant 



phase (f we have allowed for the phase factor e"**^ |T2|p^ . Obviously, the generator 
appearing on the rhs of the stochastic Schrodinger equation is non-Hermitian for 
which reason the normalization of the vector \iljg{t)) is not preserved in time. 

Due to its linearity in the state vector \iljg{t)) and the bilinearity of the generator 
in a and a^ , the stochastic Schrodinger equation ( p.l|) can be solved rigorously. For 
instance, initial coherent states | ± a), defined as eigenstates of a with eigenvalues 



±a, evolve into 21,22 



c±(t)|±ae-i*), (2.2) 

i. e. coherent states of amplitudes ±ae~2* multiplied with stochastic factors of non- 
unit modulus. 



c±(t) = exp |i («2 + |«|2) (e-^* - l)±a^j\-i'de{s)\. 



(2.3) 



Of special interest to us will be an initial superposition of two coherent states, 

mt = 0)) = ^i\a) + \-a)). (2.4) 

For large values of the amphtudes, |a| ^ 1, for which the constituents | ± a) be- 
come "macroscopically distinguishable", such a superposition is customarily called 
a Schrodinger cat state. Invoking the superposition principle we obtain the state 
originating from the superposition (|2.4|) as 



\Mt)) = ^ (c+(t)|«e-i*) + c_(t)| - ae~i*)) . (2.5) 

We can now address the mean of the quadrature component 

X^ = l (ae-^'^ + a^e'"^) (2.6) 

with respect to the normalized version of the time-evolved cat state ( p.5| ). To avoid 
unnecessary complication we shall henceforth take the amplitude a as real and thus 
get 

{x^{t))g = Arg-\Mt)\x^\Mt)) 

= Me-\a/2)e-'^4 (|c+p - |c_|') cos(^ 

"C J- T J I 

+ i{c\c^ — c+c^)(ae~2 I — a;e~2 ) sin^? >, 
Are = {Mt)\Mt))- (2.7) 

Clearly, this mean displacement vanishes initially, due to the symmetry of the ini- 
tial state, and for large times, t ^ I/7, due to the decay of lipeit)). However, at 
intermediate times we encounter a nonvanishing functional of the noise 0{t). The 
dependence of {X^p{t))g on the phase of the local oscillator is interesting: For ip = 
the overlap of the two normalized coherent states | ± ae" 2*), 



(ae ■^ \ — ae ^ 



) = exp [-2a^e-^') , (2.8) 



enters only through the normahzation factor Me'-, for ip = 7r/2, on the other hand, 
the presence of that overlap is essential for the transient appearance of a finite mean 
{XiT/2{t))e, the exponential smallness of {a\ — a) = exp(— 2«^) notwithstanding. 
Having pointed out this delicate dependence on the phase ip we proceed to setting 
(p = and writing X instead of Xq, for simplicity. 

We shall mostly be concerned with times small compared to the "mechanical" 
relaxation time I/7. In that regime the mean displacement takes the simple form 

{X{t))e = atanh{2{a^)eit)), t < I/7; (2.9) 

this results by approximating as e~ 2* ^ 1 and dropping the exponentially small over- 
lap {a\ — a). We should note that the conditional mean {X{t))e, being a functional 
of the noise ^(s) in < s < t, is a random quantity itself. It fluctuates from one 
realization of the process 6(t) to another one, each of which one might feel tempted 
to associate with a particular run of an experiment. 

In order to compute full quantum averages of oscillator observables Y like the 
displacement X we still have to average over the noise d{t) which we now need to 
specify. We want to describe an oscillator weakly interacting with a zero-temperature 
reservoir. By invoking the usual Born and Markov approximation, the stochastic 



Schrodinger equation was derived in [p!7| , pT] , p2| , together with the prescription for 
calculating quantum means 



{Yit))=Jdf^rm){Mt)\Y\Mt))- (2.10) 

We here encounter the Wiener measure (i/ij^({^}) for finding a realization 6{s) of the 
noise during the time interval < s < t. The result ( p.lC ) allows us to think of 6{t) 



as the vacuum fluctuations of the reservoir forced on the oscillator as an "input" . 

An alternative useful way of thinking of the noise is suggested by introducing the 
normalized version of the state I'lpeit)) 
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{Y{t))0 = {Mt)\Y\Mt))/{Mt)\Mt)) (2.11: 

and rewriting ( |2.1(]| ) as 



{yit))= |df^rm){Mt)\Mt)){Yit))e 

'dMm{Y{t))e. (2.12) 



We may interpret 

dfitHe}) = df,Y{{0}){Mt)\Mt)) (2.13) 

as a noise measure as well ||T5|- p^] ; indeed, by choosing F = 1 in ( |2.12|) we have the 



correct normalization J diJ,t{{6}) = 1. As was noted in |[T7|j2l|j2^ one may regard the 
random process with the measure d^t{d) as the output of the oscillator driven by 
vacuum fluctuations, i. e. a Wiener process dS,(t) as input, 

de{t) = 2^{X^{t))gdt + diit). (2.14) 

Here, we have intentionally restored the index ip on the quadrature component X^n 



since the input-ouput connection ( p.l4|) remains in fact valid for arbitrary values of 
the phase of the local oscillator. 

The two interpretations of the noise suggest different possibilities of calculating 
the means (Y{t)) by numerical simulation. One, and in fact the simpler one, is to 
generate random numbers with the Wiener measure d^Y ■> ^se them in evaluating 
lipeit)), and then employ (|2.1CI| ). According to the second strategy one would again 



construct realizations of the Wiener process and a set of ensuing |V'e(^)) but then 
proceed to the new measure d^it and finally average as required in ( p.l2| ). We should 
stress that the expectation value {X{t))g in the definition ( p. 14 ) is meant in the sense 



of the second strategy, i. e. conditioned on the measure dfif for the past history 9{s) 
s<t. 



We now propose a closer inspection of the measure d^t{d) for times < t ^ I/7. 
As already noted above we may 



then drop the overlap (ae 



(V^e(t)|V^e(t)) = -(|c- 



ae -2 ) K, l^a\ — a) , thus obtaining 



c_ 



-2a^yt 



cosh 



2a^ f de{s) 



t< 1/ 



(2.15) 



In writing out the measure dfit{{0}) niore explicitly it is convenient to discretize the 
time as tn = riAt, n = . . . N, t^ = t and to consider increments A„6' such that 
diU) = Er=o^ AiO. Assuming ^(0) = we then have 



dMO}) = \ 



'■ exp (- (A,0 - 2a^/\tf /2At 



IX V2nAt 

dA,e 



N-l 



+ n \ exp (- (A,^ + 2a^/\tf /2At 



(2.16) 



When averaging the mean displacement ( p.9|) with this measure we may first reduce 
d^t{{d}) to the marginal density of the single variable 9{t) = 9 



Pie,t) = Jdfiti{e})6(e-Y.A,e\ 



p(±) 



i{pW(e,t) + p(-)(0,t)}, 



exp{-(^T2(aV7)i)V2t}. 



^27ri 



(2.17) 



The two Gaussians here appearing are not resolved from one another as long as their 
separation 4:{ay/j)t is smaller than their width \/i. However, as soon as 

1 



i> 



27^2 



t, 



dec 



(2.18) 



the density P{6, t) has two non-overlapping peaks. Since the so-called decoherence 
time 1/270^ p7| , |28[| is, for strongly excited coherent states (a^ ^ 1), much shorter 
than the mechanical life time of excitations I/7, the separation in question takes 
place well within the range of validity (t ^ I/7) of the various simplifications made 
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in arriving at ( |2.17|) . We may conclude that the displacement of the oscillator has 
an overwhelming likelihood, at all times in the interval 1/270;^ ^ t ^ I/7, to take 
on either one or the other of two "macroscopically" distinguishable values rather 
than their vanishing ensemble average. Fig. |1| shows the density P{6,t) for various 
moments of time. 

Moreover, one intuitively expects jumps between the two macroscopically distinct 
values of {X{t))0 to become less and less frequent as t begins to exceed the deco- 
herence time tdec- In order to substantiate this expectation we consider the joint 
density 

7V-1 \ / 7V'-1 



p{e,t;e',t') = Jdfit,{{e})6l0 - ^ a A sh'-Yl ^jd 



(2.19) 



with t' = N'At > t = NAt. By noting 9' - 9 = EZn^ ^iO and doing the Gaussian 
integrals we arrive at 



P{9,t;9',t') = -lp^+\9,t)P^+\9' -9,t' -t) 

+P^-\9,t)P^~\9'-9,t'-t) 



(2.20) 



It is interesting to realize how the double-peak structure ( |2.17| ) of the single- 
point density generalizes to the two-point density (|2.2CI|) . In particular, the absence 
of cross terms P^+^P^^^ suggests that there is no interference between what happens 
near {X)g ^ +a and near (X)^ ^ —a. For a quantitative estimate of the frequency 
of jumps between the two branches we measure t and t' in units of the decoherence 
time (t = r/27Q;^) and evaluate the probability of finding 9 within a width y^ around 
the peak at +2(0:^/7)^, and 9' within a width y/f around the peak at — 2(a^/7)t', 



p{t, t') = prob 

and 9' 



,{,( 



r 



27^2^ 
/ t' 



27^2 






., (2,21) 



We immediately obtain this "jump probability" as 

p{t,t') = — lf dxe-^' I dye-y" + f rfxe"^' / dy 6"^'] (2.22) 

27r [Jii JI2 Ji[ Ji'^ J 

with the integration intervals 

Ii = [-1/V8, +1/V8] 



I 



\/t — T 

l\= [2V7-l/y8,2v^+l 



-2r' -Tx- Jt'/8 , -2r - tx + Jt'/8 



I'o 



-TX - \/t'/8 , -rx + yrV^ 



(2.23) 



Vr' — r 

Now first consider r ~ 1, r' ^ 1. Then the first summand in the jump probability 
( p.22| ) becomes negligibly small since I2 covers a relatively small interval far out in the 



wing of the Gaussian integrand; the second summand in ( |2.22| ), however, approaches a 
limit independent of r' since I2 -^ [— 1/v^, +l/v^]; therefore, the jump probability 
p{t ^ l,r' ^ 1) is not small compared to unity. However, now consider what 
happens for large times, r' > r ^ 1; then both summands in ( p.22| ) tend to become 
vanishingly small, i. e. of the order e~^^ as r grows large. Indeed, once one has 
waited several decoherence times after the preparation of the cat state, jumps between 
the macroscopically distinguishable values {X)0 = ±a become exceedingly unlikely. 
This behavior is fully borne out by numerical simulations, as exemplified in Fig. |^ 
where two "trajectories" {X{t))g pertaining to two realizations of the noise 6{t) are 
displayed. 

The foregoing reasoning makes for a certain temptation to associate single re- 
alizations of the noise 6{t) and the accompanying {X{t))g with single runs of an 
experiment. We shall discuss the legitimacy of such an interpretation further be- 
low. It may be appropriate, however, to right away ease the temptation mentioned 
by throwing a glance at Fig. ^ which depicts a simulation of (X^/2(i^))e for a sin- 
gle realization of the noise 0{t). Obviously, in contrast to (Xo(t))6) the mean of the 
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quadrature component X^/2 shows no tendency towards choosing one of the two 
distinguished values ±aexp(— ^t); even until times of the order of the mechanical 
relaxation time I/7 the mean (X^/2(i^))e stays close to zero before beginning to dis- 
play some sizable fluctuations and then eventually quieting down at zero again for 
t ^ 1/7. This behavior pO| , |3T|] of {X.„i2{t))e is easily understood by recalling the 



remark on the dependence of that mean on the phase (p of the local oscillator, made 
after (p.7|) ; for ip = 7r/2 one must wait until the overlap ( p.8|) of the two coherent 
states I ± aexp(— ^t)) has grown to values of order unity before one can expect to 
see sizable fluctuations of the mean {XT^/2{t))e away from zero. 

III. EVOLUTION OF CAT STATES ACCORDING TO THE NONLINEAR 

SCHRODINGER EQUATION 



In |T5|JT7| , |2T|j2^ the linear stochastic Schrodinger equation ( p.l| ) was shown to be 



equivalent to the following nonlinear stochastic Schrodinger equation 

d\Mt)) = { - |(«^« - 2a(X(t))^ + {Xit))l)dt 

+^{a-{Xit))^)dm}\Mt)) (3.1) 

where the mean displacement {X(t))^ = ^(a + a"'")^ is meant with respect to the state 
\(j)^) itself, 

(X(t))^ = (0^(t)|i(a + at)|0^(t)). (3.2) 

The noise increment d^{t) represents a Wiener process. The derivation of the non- 
linear evolution equation ( p.l|) from the linear one, ( p.l|) proceeds in two steps 
||r5| , p!7|pl| , |22[| . First, one interprets the noise d6{t) in (|2.1| ) as output noise driven 



by a Wiener input d^{t) according to (|2.14|) . By then normalizing as 

l<^?) = iMiO)) / \/ iM{0)\M{0)) (3-3) 
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one arrives at (|3.1| ). This equation has been employed previously by several authors 
(see e. g. ^MM)- 

Due to the equivalence of the linear and the nonlinear stochastic Schrodinger equa- 
tion and due to the existence of a rigorous solution of the hnear equation, it should 
also be possible to construct the exact solution of the nonlinear equation. This has 
in fact been achieved recently by Carmichael, Kochan and Tian [^] in the following 



way. One removes part of the nonlinearity from the nonlinear stochastic Schrodinger 
equation (|3.1| ) by working with the non-normalized wave vector \i'e{{c}){'t)) (see (|3.3|) ), 
thus obtaining 

d\Mmit)) = { - 1(«^« - MX{t))^dt + ^adm]\i^em)it))- (3-4) 

The ansatz 

\^em){t)) = e°'(^"''-^) (e^W|ae-2*) + e-'^Wj - ae'^*)) (3.5) 

yields an Ito differential equation for the random process x{t)- The associated Fokker- 
Planck equation for the density P(x, t) happens to allow for the solution given in 
i30| . After simplifying as before, i. e. e~'^* — 1 ^ — 7t, {a\ — a) ^ 0, the solution 
P(x, t) originating from P(x, 0) = 6{x) takes the form of a sum of two Gaussians, 



which in structure resembles our density (|2.17|) of the output noise 9{t). It is most 



interesting to again see the width of each of the two Gaussians grow as \/t while 

their separation grows as t, i. e. much faster. Moreover, we again find the separation 

to become manifest after a time of the order of the decoherence time tdec = l/27a^. 

The equivalence of the two versions (O) and ( |3.1| ) of the stochastic Schrodinger 



equation is worth one more illustration. The nonlinear equation ( |3.1| ) is easily seen 
to imply the following equation for the mean (X(t))j defined in 



12 



d{X{t))^ = 2^{o? - {X{t))l)dm- (3.7) 

The Ito integral of this equation takes the form 

(X(t))g = atanh (2{a^)^{t) + 47a f dt' {X{t'))^ . (3.8) 

Obviously, we have arrived at an integral equation rather than an explicit solution. 
However, the equivalence of the integral equation (|3.8|) with the solution (^]9|) of the 
linear Schrodinger equation can clearly be seen with the help of eq. (|2.14|) . 

IV. PREPARATION OF A CAT STATE 

Imagine a spin 1/2 (or equivalently, a two- level atom) prepared in an eigenstate 
of the component Sx with eigenvalue +h/2. In the Sz representation that state will 
take the form 

72(1+) + l~) ) "with Sz\±) = ±||±). Let, on the other hand, an harmonic 
oscillator be prepared in the vacuum state |0), with a\0) = 0. If we couple the two 
systems impulsively according to the Hamiltonian 

H{t) = 6{t)i2{aa^ - a*a)Sz (4.1) 

we will produce the composite state 

1 



^ya)\+) + \-a)\-)) (4.2) 

in which the coherent state |a) of the oscillator is correlated with the "up" state |+) 
of the spin while the coherent state of the opposite amplitude, | —a), is correlated 
with the spin-down-state |— ). Needless to say, the composite state ( |4.2| ) is not a cat 
state of the oscillator, even if the amplitude a is large. 

In order to proceed towards preparing a cat state we let the spin be exposed to a 
magnetic field in the y direction such that spin states are transformed as 
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l±)-e-^-«/2^|±) = -^(|+)±|-)). (4.3) 

This corresponds to a spin rotation by 7r/2 about the y axis. The composite state 
( [4.2|) is thus transformed into 

;^{W;^(|a> + |-«)) + |-);^(|a>-|-«))). (4.4) 

Now both eigenstates of 5*2 are correlated with cat states of the oscillator, if |a| ^ 1. 
We may finally imagine a measurement of Sz- Every time we find 5*^ = +h/2 we 
know that the oscillator is prepared in the cat state 

i(|«) + !-«)). 

In principle, the thought experiment just sketched can be realized by playing with 
two level atoms traversing microwave or optical cavities. Care would have to be taken 
that the two-step preparation takes less time than a decoherence time tdec = l/27|ap 
where 7 is the damping constant of the cavity within which the cat state is to be 
produced. 

V. CONCLUSION AND DISCUSSION 

Mostly employing a linear version of the stochastic Schrodinger equation we have 
presented a rigorous description of the fate of an initial state of the cat type for the 
damped harmonic oscillator. While our description is stochastically equivalent to 
one using the master equation for the density operator, it has the additional appeal 
of being more closely related to individual runs of an experiment rather than to 
ensembles of such. Indeed, the fate of the cat state described here makes for a certain 
temptation to associate a single realization of the noise 9{t) and the accompanying 
mean displacement {X{t))g with a single run of an experiment. In fact, such an 
interpretation is known to be legitimate only with respect to balanced homodyne 
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experiments where the output of the quantum oscillator is put to interference with 
a local oscillator of large amplitude and stable phase ip = 0; sharp values are then in 
effect recorded for the output noise 9(t) of the quantum oscillator. 

Inasmuch as the mean displacement {X{t))g tends to assume, for times in the 
interval 1/70^ ^ t ^ I/7, either one of the two values ±a, one may even feel 
inclined towards regarding the displacement X itself as a measured observable. The 
exceedingly small likelihood of jumps between the two preferred values ±a during 
the time interval mentioned does permit some such indulgence. One must keep in 
mind, though, that both for early times, t ^ 1/70^, and large times, t ~ I/7, the 
displacement cannot be considered as sharp in an individual run of the experiment; 
nor can the quadrature component XT^/2{t) at any time. 

Some further caution is indicated against hurried conclusions for quantum mea- 
surement theory. The heat bath providing the damping of the quantum oscillator 
is, strictly speaking, not the one needed for "objectivation" of sharp "pointer read- 
ings"; nor is the displacement X really playing the role of a pointer variable. In 
our context the quantum oscillator rather acts as a measured object and objectiva- 
tion only arises in the photodetector which generates a macroscopic electric current. 
The pointer variable with respect to which the detector must secure decoherence, 
i. e. objectivation, is just the output noise 9(t) for which sharp values are recorded. 
A striking manifestation of the potential "nonobjectivity" of the interaction of the 
oscillator with the heat bath may also be seen in the sensitive dependence of the 
statistics of X^(t) on the phase (p of the local oscillator, discussed in Section and 
already noted in ||30|j31[l . Astounding as that sensitivity may be, it does not make for 



worries in the context of quantum measurement theory; while the phase of the local 
oscillator does not influence the interaction of the quantum oscillator with its damp- 
ing reservoir, that phase does enter the oscillator observable coupled to the detector. 
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i. e. the ouput noise 9(t). The interested reader is referred to Zurek's discussion of 



"pointer bases" with respect to which fast decoherence takes place [25 
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FIGURES 
FIG. 1. Marginal density of the output noise 9 for different times r measured in units 

of the decoherence time for parameter values 7 = 0.1, |q| = 3. 

FIG. 2. The mean {XQ{t))g for two runs of the nonlinear stochastic Schrodinger equa- 
tion with initial amplitudes a = ±3. The dashed lines show the mechanical relax- 
ation iba exp(— 7t/2) of the amplitudes of the coherent states. In (a), the approach to 
+a;exp(— 7t/2) at t « 0.75 tdec is not quite close enough to produce definite locking which 
occurs only at i f« 2.6 tdcc- 

FIG. 3. The mean {X^/2{t))g. Again, a = ±3. 
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